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_, Abstract 

i-^ \ We present a pedagogical introduction to self-organized criticality (SOC), 

C ' unraveling its connections with nonequilibrium phase transitions. There are 

fi ' several paths from a conventional critical point to SOC. They begin with an 

absorbing-state phase transition (directed percolation is a familiar example), 
CN , and impose supervision or driving on the system; two commonly used methods 

^ \ are extremal dynamics, and driving at a rate approaching zero. We illustrate 

l/^ ' this in sandpiles, where SOC is a consequence of slow driving in a system 

^^ ■ exhibiting an absorbing-state phase transition with a conserved density. Other 

paths to SOC, in driven interfaces, the Bak-Sneppen model, and self-organized 
0^ , directed percolation, are also examined. We review the status of experimental 

realizations of SOC in light of these observations. 
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I. INTRODUCTION 

The label "self-organized" is applied indiscriminately in the current literature to order- 
ing or pattern formation amongst many interacting units. Implicit is the notion that the 
phenomenon of interest, be it scale invariance, cooperation, or supra-molecular organization 
(e.g., micelles), appears spontaneously. That, of course, is just how the magnetization ap- 
pears in the Ising model; but we don't speak of "self-organized magnetization." After nearly 
a century of study, we've come to expect the spins to organize; the zero-field magnetization 
below Tc is no longer a surprise. More generally, spontaneous organization of interacting 
units is precisely what we seek, to explain the emergence of order in nature. We can expect 
many more surprises in the quest to discover what kinds of order a given set of interactions 
lead to. All will be se//-organized, there being no outside agent on hand to impose order! 

"Self-organized criticality" (SOC) carries greater specificity, because criticality usually 
does not happen spontaneously: various parameters have to be tuned to reach the critical 
point. Scale- invariance in natural systems, far from equilibrium, isn't explained merely 
by showing that the interacting units can exhibit scale invariance at a point in parameter 
space; one has to show how the system is maintained (or maintains itself) at the critical 
point. (Alternatively one can try to show that there is generic scale invariance, that is, that 
criticality appears over a region of parameter space with nonzero measure |I|,0].) "SOC" has 
been used to describe spontaneous scale invariance in general; this would seem to embrace 
random walks, as well as fractal growth ||^, diffusive annihilation {A + A ^ and related 
processes), and nonequilibrium surface dynamics |^. Here we restrict the term to systems 
that are attracted to a critical (scale-invariant) stationary state; the chief examples are 
sandpile models [^. Another class of realizations, exemplified by the Bak-Sneppen model 
0], involve extremal dynamics (the unit with the extreme value of a certain variable is the 
next to change). We will see that in many examples of SOC, there is a choice between 
global supervision (an odd state of affairs for a "self-organized" system), or a strictly local 
dynamics in which the rate of one or more processes must be tuned to zero. 

The sandpile models introduced by Bak, Tang and Wiesenfeld (BTW) 0], Manna 0, 
and others have attracted great interest, as the first and clearest examples of self-organized 
criticality. In these models, grains of "sand" are injected into the system and are lost at the 
boundaries, allowing the system to reach a stationary state with a balance between input 
and output. The input and loss processes are linked in a special way to the local dynamics, 
which consists of activated, conservative, redistribution of sand. In the limit of infinitely 
slow input, the system displays a highly fluctuating, scale-invariant avalanche-like pattern 
of activity. One may associate rates h and e, respectively, with the addition and removal 
processes. We have to adjust these parameters to realize SOC: it appears in the limit of h 
and e — i> 0"*" with h/e — >■ [pl,PHTO[]. (The addition and removal processes occur infinitely 
slowly compared to the local redistribution dynamics, which proceeds at a rate of unity. 
Loss is typically restricted to the boundaries, so that e — *> is implicit in the infinite-size 
limit.) 

Questions about SOC fall into two categories. First, Why does self-organized criticality 
exist? What are the conditions for a model to have SOC? Second, the many questions about 
critical behavior (exponents, scaling functions, power-spectra, etc.) of specific models, and 
whether these can be grouped into universality classes, as for conventional phase transitions 



both in and out of equilibrium. Answers to the second type of question come from exact 



solutions 1 11], simulations [12|, renormalization group analyses |13|, and (one may hope) 
field theoretical analysis. Despite these insights, assertions in the literature about sponta- 
neous or parameter-free criticality have tended to obscure the nature of the phase transition 
in sandpiles, fostering the impression that SOC is a phenomenon sui generis, inhabiting 
a different world than that of standard critical phenomena. In this paper we show that 
SOC is a phase transition to an absorbing state, a kind of criticality that has been well 
studied, principally in the guise of directed percolation |T^. Connections between SOC and 



an underlying conventional phase transition have also been pointed out by Narayan and 
Middleton [|T^, and by Sornette, Johansen and Dornic |I6| . 

Starting with a simple example (Sec. II), we will see that the absorbing-state transition 
provides the mechanism for SOC (Sec. III). That is, we explain the existence of SOC in 
sandpiles on the basis of a conventional critical point. In Sec. IV we discuss the trans- 
formation of a conventional phase transition to SOC in the contexts of driven interfaces, 
a stochastic process that reproduces the stationary properties of directed percolation, and 
the Bak-Sneppen model. We find that criticality requires tuning, or equivalently, an infinite 
time-scale separation. With this essential point in mind, we present a brief review of the rel- 
evance of SOC models to experiments in Sec. V. Sec. VI presents a summary of our ideas. 
We note that this paper is not intended as a complete review of SOC; many interesting 
aspects of the field are not discussed. 

II. A SIMPLE EXAMPLE 

We begin with a simple model of activated random walkers (ARW). Each site j of a lattice 
(with periodic boundary conditions) harbors a number Zj = 0, 1, 2... of random walkers. (For 
purposes of illustration the ring 1, ..., L will do.) Initially, A^ walkers are distributed randomly 
amongst the sites. Each walker moves independently, without bias, to one of the neighboring 
sites (i.e., from site j to j ' + 1 or j — 1, with site L + 1 = 1 and = L), the only restriction 
being that an isolated walker (at a site with Zj = 1) is paralyzed until such time as another 
walker or walkers joins it. The active sites (with Zj > 2) follow a Markovian (sequential) 
dynamics: each active site loses, at a rate 1, a pair of walkers, which jump independently 
to one of the neighbors of site j. (Thus in one dimension there is a probability of 1/2 that 
each neighbor gains one walker, while with probability 1/4 both walkers hop to the left, or 
to the right.) 

The model we have just defined is characterized by the number of lattice sites, L'^, and 
the number of particles, A^. It has two kinds of configurations: active, in which at least one 
site has two or more walkers, and absorbing, in which no site is multiply occupied, rendering 
all the walkers immobile [17|. For N > L'^ only active configurations are possible, and since 



A^ is conserved, activity continues forever. For N < L'^ there are both active and absorbing 
configurations, the latter representing a shrinking fraction of configuration space as the 
density C = N/ L'^ —> 1. Given that we start in an active configuration (a virtual certainty 
for an initially random distribution with C > and L large), will the system remain active 
indefinitely, or will it fall into an absorbing configuration? For small ( it should be easy for 
the latter to occur, but it seems reasonable that for sufficiently large densities (still < 1), the 
likelihood of reaching an absorbing configuration becomes so small that the walkers remain 



active indefinitely. In other words, we expect sustained activity for densities greater tfian 
some critical value (c, with (^ < ^■ 

A simple mean-field theory provides a preliminary check of this intuition. Consider 
activated random walkers in one dimension. For a site to gain particles, it must have an active 
{z > 2) nearest neighbor. Since active sites release a pair of walkers at a rate of unity, a given 
site receives a single walker from an active neighbor at rate 1/2, and a pair of walkers at rate 
1/4. Thus the rate of transitions that take Zj to Zj + 1 is [P{zj, Zj^i > 2) + P{zj, Zj-i > 2)]/2; 
transitions from Zj to Zj + 2 occur at half this rate. In the mean-field approximation we 
ignore correlations between different sites, and factorize the joint probability into a product: 
P{z, z' > 2) = PzPa, where p^ is the fraction of sites with occupation z and pa = J2z>2 Pz is 
the fraction of active sites. Using this factorization, we can write a set of equations for the 
site densities: 

-^ = Pa{pz-1- Pz) + ^Pa{pz-2- Pz)+Pz+2-0^_2pz , {z = 0,1,2...), (l) 

where ^^ = for n < and is one otherwise. The final two terms represent active sites losing 
a pair of walkers. It is easy to see that the total probability, and the density ( = J2z ^Pz 
are conserved by the mean-field equations. This infinite set of coupled equations can be 
integrated numerically if we impose a cutoff at large z. (This is justified by the finding 
that Pz decays exponentially for large z.) The mean-field theory predicts a continuous phase 
transition at Cc = 1/2- For ( < (c the only stationary state is the absorbing one, pa = 0, 
while for ( > (c the active-site density grows oc C — Cc- A two-site approximation (in which 
we write equations for the fraction p^^z' of nearest-neighbor pairs with given heights, but 
factorize joint probabilities involving three or more sites), yields Cc = 0.75. 

The existence of a continuous phase transition is confirmed in Monte Carlo simulations, 
which yield Cc — 0.9486 in one dimension, and Cc — 0.7169 in two dimensions. Figure 1 shows 
how the stationary density of active sites pa depends on (; we see pa growing continuously 
from zero at (c- (The points represent estimated densities for L — >■ cxd, based on simulation 
data for L = 100 — 5000.) The inset shows that the active-site density follows a power law. 
Pa ~ (C ~ Cc)^5 with P = 0.43(1); a finite-size scaling analysis confirms this result [0. In 



summary, activated random walkers exhibit a continuous phase transition from an absorbing 
to an active state as the particle density is increased above Cc with (c strictly less than 1. 
(It has yet to be shown rigorously that the active-site density in the ARW model is singular 
at (c, in the infinite-size limit; our numerical results are fully consistent with the existence 
of such a singularity.) 

A. Absorbing-State Phase Transitions 

Absorbing-state phase transitions are well known in condensed matter physics, and pop- 



ulation and epidemic modeling |]T9l. The simplest example, which may be thought of as the 
"Ising model" of this class of systems, is the contact process [^ . Again we have a lattice of 
L'^ sites, each of which may be occupied (active) or vacant. Occupied sites turn vacant at 
a rate of unity; vacant sites become occupied at a rate of {X/2d)na where n^ is the number 
of occupied nearest neighbors (the factor 2d represents the number of nearest neighbors). 
There is a unique absorbing configuration: all sites vacant. For A sufficiently small, the 



system will eventually fall into the absorbing state, while for large A an active stationary 
state can be maintained. Letting p represent the density of occupied sites, the mean-field 
theory analogous to the one formulated above for activated random walkers reads: 

|^(A-1)P-A/. (2) 

This predicts a continuous phase transition (from p = Otop=l — A~^ in the stationary 
state) at Ac = 1. Rigorous analyses [^,0 confirm the existence of a continuous phase 



transition at a critical value Ac, in any dimension d > 1. Simulations and series analyses 
yield Ac = 3.29785(2) in one dimension. This model, and its continuous- up date counterpart, 
directed percolation (DP; see Sec. IV), have been studied extensively. The critical exponents 
are known to good precision for d = 1, 2, and 3; the upper critical dimension dc = 4. There 



is, in addition, a well established field theory for this class of models [23,24|: 



^ = V^p-ap-bp^+r]{x,t). (3) 

Here p{x,t) is a local particle density, and ri{x,t) is a Gaussian noise with autocorrelation 

{r]{x, t)r]{x', t')) = Tp{x, t)6{x - x')6{t - t') . (4) 

That (r/^) is linear in the local density follows from the fact that the numbers of events 
(creation and annihilation) in a given region are Poissonian random variables, so that the 
variance equals the expected value. (The noise must vanish when p = for the latter to be 
an absorbing state!) This field theory serves as the basis for a strong claim of universality 
23| , |25[| : Continuous phase transitions to an absorbing state fall generically in the universality 



class of directed percolation. (It is understood that the models for which we expect DP- 
like behavior have short-range interactions, and are not subject to special symmetries or 
conservation laws beyond the simple translation-invariance of the contact process. Models 



subject to a conservation law are known to have a different critical behavior ||26|| .) 

The activated random walkers model resembles the contact process in having an 
absorbing-state phase transition. We should note, however, two important differences be- 
tween the models. First, ARW presents an infinite number (2^ , to be more precise) of 
absorbing configurations, while the CP has but one. In fact, particle models in which the 
number of absorbing configurations grows exponentially with the system size have also been 
studied intensively. The simplest example is the pair contact process, in which both elemen- 
tary processes (creation and annihilation) require the presence of a nearest-neighbor pair of 
particles |^. In one dimension, a pair at sites i and i + 1 can either annihilate, at rate 
p, or produce a new particle at either i — 1 or i + 2, at rate 1 — p (provided the selected 
site is vacant). This model shows a continuous phase transition from an active state for 
p < Pc to an absorbing state above pc- The static critical behavior again belongs to the DP 
universality class, but the critical exponents associated with spreading of activity from an 
initially localized region are nonuniversal, varying continuously (in one dimension) with the 



particle density in the surrounding region [28 



A second important difference between ARW and the CP and PCP is that the former is 
subject to a conservation law (the number of walkers cannot change from its initial value). 
In a field-theoretic description of ARW we will therefore need (at least) two fields: the local 



density p{x,t) of active sites, and the local particle density ({x,t); the latter is frozen in 
regions where p = 0. The evolution of p is coupled to ( because the particle density controls 
existence and level of activity in the ARW model. 

Given that absorbing-state phase transitions fall generically in the universality class of 
directed percolation, it is natural to ask whether this is the case for activated random walkers 
as well. The answer, apparently, is "No." The critical exponent (3 for ARW is, as we noted 
above, 0.43, while for one-dimensional DP (3 = 0.2765 P9; the other critical exponents differ 



as well |]T8[. While the reason for this difference is not understood, it appears, at least, to 
be consistent with the existence of a conserved field in ARW. 

To summarize, our simple model of activated random walkers has an absorbing-state 
phase transition, as does the contact process, directed percolation and the POP. All possess 
the same basic phase diagram: active and inactive phases separated by a continuous phase 
transition at a critical value of a "temperature-like" parameter {( in ARW, A in the CP). 
But ARW possesses an infinite number of absorbing configurations, and the evolution of its 
order parameter (the active-site density) is coupled to a conserved density (. The latter 
presumably underlies its belonging to a different universality class than DP. 

III. ACTIVATED RANDOM WALKERS AND SANDPILES 

The activated random walkers model possesses a conventional critical point: we have to 
tune the parameter ( to its critical value. What has it got to do with self- organized criti- 
cality? The answer is that ARW has essentially the same local dynamics as a model known 
to exhibit SOC, namely, the Manna sandpile |^. In Manna's sandpile, the redistribution 
dynamics runs in parallel: at each time step, all of the sites with z > 2 simultaneously 
liberate two walkers, which jump randomly to nearest neighbor sites. This may result in 
a new set of active sites, which relax at the next time step, and so on. (Time advances 
by one unit at each lattice update, equivalent to the unit relaxation rate of an active site 
in ARW.) We defined ARW with sequential dynamics as this makes it a Markov process 
with local transitions in configuration space, like a kinetic Ising model. There is of course 
nothing wrong in defining ARW with parallel dynamics; it too has an absorbing-state phase 
transition. 

There is a much more fundamental difference between the Manna sandpile and the ARW 
model: the former allows addition and loss of walkers. Recall that we defined the ARW with 
periodic boundary conditions; walkers can never leave the system. In the sandpile walkers 
may exit from one of the boundary sites. (On the square lattice, for example, a walker at an 
edge site has a probability of 1/4 to leave the system at the next step.) If we allow walkers to 
leave, then eventually the system will reach an absorbing configuration. When this happens, 
we add a new walker at a randomly chosen site. This innocent- sounding prescription — 
add a walker when and only when all other activity ceases — carries the infinite time 
scale separation essential to the appearance of SOC in sandpiles. The sequence of active 
configurations between two successive additions is known as an avalanche; avalanches may 
involve any number of sites, from zero (no topplings) up to the entire system. 

Manna showed that his model reaches a stationary state in which avalanches occur on 
all scales, up to the size of the system, and follow a power-law distribution, P{s) ~ s~'^, 
for s -C Sc- (Here s is the number of transfer or toppling events in a given avalanche, and 
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Sc ~ L^ is a cutoff associated with the finite system size.) In other words, the Manna 
sandpile, hke the models devised by Bak, Tang and Wiesenfeld and others, exhibits scale 
invariance in the stationary state. 

We know that ARW, which has the same local dynamics as the Manna sandpile, shows 
scale invariance when (and only when) the density ( = Cc- So in the stationary state 
of the Manna model, the density is somehow attracted to its critical value. How does it 
happen? The mechanism of SOC depends upon a particular relation between the input 
and loss processes, and the conventional absorbing-state phase transition in the model with 
a fixed number of particles. Walkers cannot enter the system while it is active, though 
they may of course leave upon reaching the boundary. In the presence of activity, then, 
( > (c and d(/dt < 0. In the absence of activity there is addition, but no loss of walkers, 
so C < Cc implies d(/dt > 0. Evidently, the only possible stationary value for the density 
in the sandpile is (J Of course, it is possible to have a low level of activity locally, in a 
region with ( < (c, but under such conditions activity cannot propagate or be sustained. 
(One can similarly construct absorbing configurations with ( > (c, but these are unstable to 
addition of walkers, or the propagation of activity from outside.) In the infinite-size limit, 
the stationary activity density is zero for ( < (c, and positive for ( > Q, ensuring that ( is 
pinned at (c, when loss is contingent upon activity, and addition upon its absence. 

That the Manna sandpile, in two or three dimensions, with parallel dynamics, has a scale- 
invariant avalanche distribution is well known . Here we note that the same holds for the 
one-dimensional version, with random sequential dynamics. Figure 2 shows the probability 
distribution for the avalanche size (the total number of topplings) when we modify ARW 
to include loss of walkers at the boundaries, and addition at a randomly chosen site, when 
the system falls into an absorbing configuration. The distribution follows a power law, 
P{s) ~ s~^% over a wide range of avalanche sizes and durations; there is, as expected, an 
exponential cutoff Sc ~ L^ for events larger than a characteristic value associated with 
the finite size of the lattice. (Our best estimates are Ts = 1.10(2) and D = 2.21(1).) The 
upper inset of Fig. 2 shows that the stationary density approaches (c, the location of the 
absorbing-state phase transition, as L — »• cx). It is also interesting to note that, in contrast 
with certain deterministic one-dimensional sandpile models [^,^, the present example 



appears to exhibit finite-size scaling, as shown in the lower inset of Fig. 2. 

A. A Recipe for SOC 

The connection between activated random walkers and the Manna sandpile suggests the 
following recipe for SOC. Start with a system having a continuous absorbing-state phase 
transition at a critical value of a density (. This density should represent the global value 
of a local dynamical variable conserved by the dynamics. Add to the conservative local 
dynamics (1) a process for increasing the density in infinitesimal steps (C ~^ C + dC) when 
the local dynamics reaches an absorbing configuration, and (2) a process for decreasing the 
density at an infinitesimal rate while the system is active. Run the system until it reaches 
the stationary state; it is now ready to display scale invariance. 

Let's see how these elements operate in the Manna sandpile. We started with activated 
random walkers, which does indeed display a continuous absorbing-state transition as a 
function the density ( of walkers; this density, moreover, is conserved. To this we added the 



input of one walker (C "~^ C + V-^*^ i^ d dimensions), when the system is inactive. We then 
broke the translational symmetry of the ARW model to define boundary sites, and allowed 
walkers at the boundary to leave the system. The latter implies a loss rate d(/dt oc —L~^pf,, 
where pi, is the activity density at the boundary sites. The conditions of our recipe are 
satisfied when L ^ oo, which we needed anyway, to have a proper phase transition in the 
original model. 

Now we can examine the ingredients one by one. First, the phase transition in the original 
model should be to an absorbing state, because our input and loss steps are conditioned on 
the absence or presence of activity. Second, the temperature-like parameter controlling the 
transition should be a conserved density. So the contact process and PCP aren't suitable 
starting points for SOC, because the control parameter A isn't a dynamical variable. (To 
self-organize criticality in the CP, we'd have to change A itself, depending on the absence of 
presence of activity. But this is tuning the parameter by hand!) Third, we need to change the 
density ( in infinitesimal steps, else we will always be jumping between values above or below 
(c without actually hitting the critical density. The same thing will happen, incidentally, if 
we start out with a model that has a discontinuous transition (with attendant hysteresis) 
between an active and an absorbing state; this yields self-organized stick-slip behavior. 

The basic ingredients of our recipe are an absorbing-state phase transition, and a method 
for forcing the model to its critical point, by adding (removing) particles when the system 
is frozen (active). Following the recipe, the transformation of a conventional critical point 



to a self-organized one does not seem surprising [32 



B. Firing the Baby-Sitter 

The reader may have noted a subtle inconsistency in the above discussion. We rejected 
the contact process as a suitable candidate for SOC because changing the parameter A on 
the basis of the current state (active or frozen) amounts to tuning. Cannot the same be said 
for adding walkers in the Manna sandpile? Somehow, a dynamics of walkers entering and 
leaving the system seems more "natural" than wholesale fiddling with a parameter. But 
who is going to watch for activity, to know when to add a particle? A system managed by 
a supervisor can hardly be called "self-organized!" If we want to avoid building a supervisor 
or baby-sitter into the model, we had better say that addition goes on continuously, at rate 
h, and that SOC is realized in the limit h ^ 0^ li;lC|. (The original sandpile definitions 



have a baby-sitter. Simulations, in particular, have a live-in baby-sitter to decide the next 
move. Addition at rate /i — *■ O"*" is a supervisor-free interpretation of the dynamics [0.) 
In the recipe for SOC without baby-sitters, we replace addition (1) above with (!'): allow 
addition at rate h, independent of the state of the system, and take h -^ 0^. (There is no 
problem with the removal step: dissipation is associated with activity, which is local.) We 
pay a price when we fire the baby-sitter: there is now a parameter h in the model, which has 
to be tuned to zero. Evidently, sandpiles don't exhibit generic scale invariance, but rather, 
scale invariance at a point in parameter space. This is consistent with Grinstein's definition 
of SOC, which requires an infinite separation of time scales from the outset f^. 



C. Variations 

In certain respects, our recipe allows greater freedom than was explored in the initial 
sandpile models. There is no special reason, for example, why loss of walkers has to occur at 
the boundaries. We simply require that activity be attended by dissipation at an infinitesimal 
rate. SOC has, indeed, been demonstrated in translation-invariant models with a uniform 
dissipation rate ep when e ^ 0''" [^,0]. In the original sandpile models, addition takes place 
with equal probability at any site, but restricting addition to a subset of the lattice will still 
yield SOC. 

Our recipe allows a tremendous amount of freedom for the starting model; the only 
restriction is that it possess an absorbing-state critical point as a function of a conserved 
density. The dynamical variables can be continuous or discrete. The hopping process does 
not have to be symmetric, as in ARW. (In fact, directed hopping yields an exactly-soluble 
sandpile [^.) The model need not be defined on a regular lattice; any structure with a 
well defined infinite-size limit should do. The dynamics, moreover, can be deterministic. 
Consider a variant of the ARW model (on a rf-dimensional cubic lattice) in which a site 
is active if it has z > 2d walkers. At each lattice update (performed here with parallel 
dynamics), every active site 'topples, transferring a single walker to each of the 2d nearest- 
neighbor sites. In this case the only randomness resides in the initial configuration. But the 
model again exhibits a continuous absorbing-state phase transition as we tune the number of 
walkers per site, (. Starting with this deterministic model, our recipe yields the celebrated 
Bak-Tang-Wiesenfeld sandpile. 

As a further variation, we can even relax the condition that the order parameter is 
coupled to a conserved field [Q. The price is the introduction of an additional driving 
rate. This situation is exemplified by the forest-fire model p7| , |38| . The model is defined 



on a lattice in which each site can be in one of three states: empty, or occupied by a tree, 
either live or burning. Burning trees turn into empty sites, and set fire to the trees at 
nearest-neighbor sites, at a rate of unity. It is easy to recognize that burning trees are the 
active sites: any configuration without them is absorbing. In an infinite system, there will 
be a critical tree density that separates a phase in which fires spread indefinitely from an 
absorbing phase with no burning trees. In a finite system we can study this critical point 
by fixing the density of trees at its critical value pSf j. 

So far we have no process for growing new trees. The forest-fire propagates like an 
epidemic with immunity: a site can only be active once, and there is no proper steady state 
^ . As in sandpiles, to obtain a SOC state we must introduce an external driving field / that 



introduces a small probability for each tree to catch fire spontaneously. This driving field 
allows the system to jump between absorbing configurations through the spreading of fires. 
The latter, however, are completely dissipative, i.e., the number of trees is not conserved. 
Thus, if we want to reach a stationary state we must introduce a second external driving 
field p that causes new trees to appear. (Empty sites become occupied by a living tree at 
rate p.) In this case criticality is reached by the double slow driving condition /, p — ^ and 
//j9 — > 0. In practice, this slow driving condition is achieved by the usual supervisor, that 
stops fire ignition and tree growth during active intervals. 



D. Fixed-Energy Sandpiles 

If someone hands us a sandpile displaying SOC, we can identify the initial model in our 
recipe; it has the same local dynamics as the SOC sandpile. Thinking of the conserved ( 
as an energy density, we call the starting model a fixed-energy sandpile (FES). Thus the 
activated random walkers model introduced in Sec. II is the fixed-energy Manna sandpile, 
and the variant described in the preceding subsection is the BTW FES. Now the essential 
feature of the fixed-energy sandpile is an absorbing-state phase transition. SOC appears 
when we rig up the addition and removal processes to drive the local FES dynamics to Cc- To 
understand the details of SOC, then, we ought to try to understand the conventional phase 
transition in the corresponding fixed-energy sandpile. This is our program for addressing 
the second class of questions (about critical exponents and universality classes) mentioned 
in the Introduction. Since fixed-energy sandpiles have a simple dynamics (Markovian or 
deterministic) without loss or addition, and are translation-invariant (when defined on a 
regular lattice), they should be easier to study than their SOC counterparts. The relation 
to absorbing-state phase transitions leads to a proper identification of the order parameter 
M, and suggests a strategy for constructing a field theory of sandpiles pi|. Spreading 



exponents, conventionally measured in absorbing-state phase transitions, are related through 
scaling laws to avalanche exponents, usually measured in slowly driven systems 

IV. OTHER PATHS TO SOC 

A. Driven Interfaces 

In this section we illustrate the central idea of the preceding section — the transformation 
of a conventional phase transition to a self-organized one — in a different, though related, 
context. We begin with a single point mass undergoing driven, dissipative motion in one 
dimension. Its position H{t) follows the equation of motion 

, ^d H dH „ „ , , , s 

where M is the mass, 'jH represents viscous dissipation, F is the applied force, and Fp{H) 
is a position-dependent pinning force. In many cases of interest (i.e., domain walls or flux- 
lines) the motion is overdamped and we may safely set M = 0. The pinning force has mean 
zero {{Fp{h)) = 0) and its autocorrelation {Fp{h)Fp{h + y)) = A{\y\) decays rapidly with 
\y\; the statistical properties of Fp are independent of H. Assuming, as is reasonable, that 
Fp is bounded {Fp < Fm), we expect the motion to continue if the driving force F exceeds 
Fm- Otherwise the particle gets stuck somewhere. 

Now consider an elastic interface (or a flux line) subject to an external force, viscous 
damping, and a pinning force associated with irregularities in the surrounding medium. If 
we discretize our interface, using ifj(t) to represent the position, along the direction of the 
driving force, of the i-th segment , the equation of motion is 



dHj 

^ dt 



H,+i + H,^i - 2H,{t) +F- Fp,,{Hi), (6) 
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where the Fp^i{Hi) are a set of independent pinning forces with statistical properties as 
above. This driven interface model has a depinning transition at a critical value, F^., of the 
driving force ^^ . (Eq. (|^) describes a linear driven interface, so-called because it lacks the 
nonlinear term oc (V/i)^, familiar from the KPZ equation p|,^.) For F < F^. the motion is 
eventually arrested [dHi/dt = for all i), while for F > Fc movement continues indefinitely. 
Close to Fc there are avalanche-like bursts of movement on all scales, interspersed with 
intervals of near-standstill. The correlation length and relaxation time diverge at Fc, as in 
the other examples of absorbing-state phase transitions we've discussed above. We may take 
the order parameter for this transition as the mean velocity, v = (dHi/dt). 

To reach the absorbing-state phase transition in the driven interface model we need to 
adjust the applied force F to its critical value Fc. Can we modify this system so that it will 
be attracted to the critical state? Note that F is not a dynamical variable, any more than 
is A, in the contact process. Our sandpile recipe doesn't seem to apply here. The crucial 
observation is that we may change the nature of the driving, replacing the constant force 
F with a constraint of fixed velocity, dHi/dt = v. A finite v corresponds to a state in the 
active phase: the mean driving force (Fi)^ > Fc for f > 0. When we allow v to tend to zero 
from above, we approach the depinning transition. This limit can be attained through an 
extremal dynamics in which we advance, at a given step, only the element subject to the 
smallest pinning force |]46|j47[| . (Notice that in extremal dynamics we are directly adjusting 
the order parameter [p!6| .) 

To avoid the global supervision implicit in extremal dynamics we may attach each element 
of the interface to a spring, and move the other end of each spring at speed V. Now the 
equations of motion read 

7^ = H,+i + /J,_i - 2H,{t) + k{Vt - Hi) - Fp,,{H), (7) 

where k is the spring constant. For high applied velocities, the interface will in general 
move smoothly, with velocity H = V, while for low V stick-slip motion is likely. In the 
overdamped regime, the amplitudes of the slips are controlled by V and k, and the statistics 
of the potential. In the limit ^ ^ 0, the interface motion exhibits scale invariance; V plays 
a role analogous to h in the sandpile. (The limits V —>■ and k —>■ have a particular 
significance, since the block can explore the pinning- force landscape quasistatically.) The 
fine tuning of F to Fc in the constant-force driving has been replaced by fine tuning V to zero. 
This parameter tuning corresponds, once again, to an infinite time-scale separation. Finally, 
we note that restoring inertia (M > 0) results in a discontinuous depinning transition with 
hysteresis, resulting in stick-slip motion of the sort associated with friction [^ . 

Once again, we have transformed an absorbing-state phase transition {F = Fc) into 
SOC by driving the system at a rate approaching zero {V —>■ 0). But there appear to be 
fundamental differences between sandpiles and driven interfaces. In the sandpile, but not in 
the driven interface, the order parameter is coupled to a conserved density. The sandpile, 
moreover, does not involve a quenched random field as does the driven interface. Despite 
these apparent differences, close connections have been suggested between the two kinds of 



model [15,49-51|. We review this correspondence in the next subsection, following Ref. [51 . 
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B. Sandpiles and Driven Interfaces 

Consider the BTW fixed-energy sandpile in two dimensions; let ifj(t) be the number 
of times site i has toppled since time zero. To write a dynamics for Hi, we observe that 
the occupation Zi(t) of site i differs from its initial value, ^i(O), due to the inflow and the 
outflow of particles at this site. The outflow is given by 4ifj(t), since each toppling expels 
four particles. The inflow can be expressed as Y^nn Hj(t): site i gains a particle each time 
one of its nearest neighbors topples. Summing the above contributions we obtain: 

z,{t) = z,{0)+ Y. H,{t)-4H,{t) 

jNNi 

= z,{0) + VlH,{t), (8) 

where V|) stands for the discretized Laplacian. Since sites with Zi(t) > 4 topple at unit 
rate, the dynamics of H^ is given by 

^ = e[zm + ^lH^{t)-3] 

= Q[VlH,{t) + F - F,^i], (9) 

where dHi/dt is shorthand for the rate at which the integer- valued variable Hi{t) jumps to 
Hi{t) + 1, and G(x) = 1 for x > and is zero otherwise. In the second line, F = ( — 3 and 
Fpj = Zi{0) — (. (Recall that ( = {zi{t)) for all t.) Thinking oi Hi{t) as a discretized interface 
height, Eq. @ represents an overdamped, driven interface in the presence of columnarnoise, 
Fp^i, which takes independent values at each site, but does not depend upon Hi, as it does 
in the interface model discussed in the preceding subsection. We see from this equation that 
tuning ( to its critical value Cc is analogous to tuning the driving force to Fc. If we replace 
the discrete height Hi in Eq. (^ with a continuous field, H{x,t) (and similarly for Fp), and 
replace the B-function by its argument, we obtain the Edwards- Wilkinson surface-growth 
model with columnar disorder, which has been studied extensively |5^. The similarity 
between the present height representation and the dynamics of a driven interface suggests 
that the critical point of the BTW fixed-energy sandpile belongs to the universality class 
of linear interface depinning with columnar noise, if the rather violent nonlinearity of the 
©-function is irrelevant. (The latter remains an open question. A height representation 
for the Manna sandpile is also possible, but is complicated by the stochastic nature of the 
dynamics.) 

Applying the recipe of Sec. Ill to the driven interface, we would impose open boundaries, 
which drag behind the interior as they have fewer neighbors pulling on them; eventually the 
interface gets stuck. When this happens, we ratchet up the "force" at a randomly chosen 
site (in effect, Fpj -^ Fpj — 1 at the chosen site). The dynamics is then attracted to the 
critical point. Once again, we may trade supervision (checking if the interface is stuck) for 
a constant drive {F —>■ F + ht) in the limit h —>■ 0. 

C. Self-Organized Directed Percolation and the Bak-Sneppen Model 

Take the square lattice and rotate it by 45°, so that each site has two nearest neighbors 
in the row above, and two below. The sites exist in one of two states, "wet" and "dry." 
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The states of the sites in the zeroth (top) row can be assigned at will; this defines the initial 
condition. A site in row i > 1 is obliged to be dry if both its neighbors in row i — 1 are 
dry; otherwise, it is wet with probability p, and dry with probability 1 — p. This stochastic 
cellular automaton is called site directed percolation. Like the contact process, it possesses an 
absorbing state: all sites dry in row k implies all dry in all subsequent rows. The dynamics 
of site DP can be expressed in a compact form if we define the site variable x* to be zero 
(one) if site j in row i is wet (dry). The variables in the next row are given by 

x;+^ = e[max{r/}, min{x}_i, x).^^}} - p] , (10) 

where the r]j are independent random variables, uniform on [0,1]. If both neighbors in the 
preceding row are in state 1, x^"*"^ must also equal 1; otherwise x*"^^ = with probability p. 
Thinking of the rows as time slices, we see that site DP is a parallel-update version of the 
contact process: increasing p renders the survival and propagation of the wet state more 
probable, and is analogous to increasing A in the CP. Just as the CP has a phase transition 
at Ac, site DP has a transition from the absorbing to the active phase at Pc — 0.7054. 

We've already dismissed the contact process (and by extension DP) as starting models 
for realizing SOC via the recipe of Sec. III. Remarkably, however, it is possible to define a 
parameter-free stochastic process whose stationary state reproduces the properties of crit- 
ical DP [^^-|55|. This process, self-organized directed percolation (SODP), is obtained by 



replacing the discrete variables in Eq. ([10|) by real variables which store the value of one of 
the previous ?7*. In place of Eq. (jT^) we have simply 

x*"^^ = max{r7*,min{x*„i,x*+i}} , (11) 

Notice that parameter p has disappeared, along with the G function. Starting from a 
distribution with x^ < 1 for at least one site (but otherwise arbitrary), this process eventually 
reaches a stationary state, characterized by the probability density fi{x). One finds that fi{x) 
is zero for x < Pc (the critical value of site DP), jumps to a nonzero value (infinity, in the 
thermodynamic limit), at Pc, and decreases smoothly with x for x > Pc- The process has 
discovered the critical value of site directed percolation! 

Hansen and Roux explained how this works [^: for any p G [0, 1] the probability that 
x'j < p is p ii either or both of the neighbors in the previous time slice have values less that 
p (i.e., if the smaller of x*~\ and x*I^\ is < p), and is zero if x*~\ and x*I,l\ both exceed p. 
This is exactly how the "wet" state propagates in site DP, with parameter p, if we equate 
the events 'site j in row i is wet' and 'x* < p.' It follows that in the stationary state, 

Pt[x] <p]= iJ,{x)dx, (12) 

•^ Jo 

equals the probability P{p) that a randomly chosen site is wet, in the stationary state of site 
DP with parameter p. This explains why /i(x) = for x < Pc, and why ^{pc) is infinite in 
the infinite-size limit {dP/dp is infinite at pc) ■ The spatio-temporal distribution of DP is also 
reproduced; for example, the joint probability Pr[x* < pc, x\ < pc] decays as a power law 
for large separations \j — k\. The process effectively studies all values of p at once, greatly 
improving efficiency in simulations. Stochastic processes corresponding to other models (DP 
on other lattices, bond instead of site DP, epidemic processes) have also been devised |]5^ , |56| . 
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It seems unlikely, on the other hand, that such a real-valued stochastic process exists for 
activated random walkers or other fixed-energy sandpiles. (Of course, such a process would 
be of great help in studying sandpiles!) 

SODP doesn't fit into the same scheme as sandpiles or driven interfaces. It is a real- 
valued stochastic process that generates, by construction, the probability distribution of DP 
for all parameter values, including pc- The process itself does not have a phase transition; all 
sites are active (except those inside a sequence of I's — a configuration that will never arise 
spontaneously), since there is a finite probability for x* to change. SODP is self-organized in 
the sense that its stationary probability density has a critical singularity, without the need to 
adjust parameters. If we choose to regard SODP as an instance of SOC, we must recognize 
that the path in this case is very different from that in sandpiles or driven interfaces; the 
system is not being forced to its critical point by external supervision or driving. Rather, 
SODP is directed percolation implemented in a different (parameter-free) way. Furthermore, 
the dynamics embodied in Eq. (|TT|) seems a much less realistic description of a physical 
system than is driven-interface motion, or even the rather artificial dynamics of a sandpile 
model. In the rather unlikely event that SODP were realized in a natural system, it would 
not immediately yield a scale-invariant "signal" such as avalanches or fractal patterns. The 
latter would require a second process (or an observer) capable of making fine distinctions 
among values of x in the neighborhood of Pc- So the kind of SOC represented by SODP 
does not appear a likely explanation of scale invariance in nature. 

A (fanciful) interpretation of Eq. (p!TD is that x*- represents the "fitness" of an individual, 
which mates with its neighbor to produce an offspring that inherits the fitness of the less-fit 
parent. This offspring survives if her fitness exceeds that of an interloper, whose fitness is 
random. (It is, to put it crudely, as if an established population were constantly challenged 
by a fiux of outsiders.) Seen in this light, SODP bears some resemblance to the evolutionary 
dynamics represented, again in very abstract form, in the Bak-Sneppen model 0. Here, 
the globally minimum fitness variable, along with its nearest neighbors, is replaced by a [0,1] 
random number at each time step. (If the x* are associated with different species, then the 
appearance of a new species at site i affects the fitness of the "neighboring" species in the 
community in an unpredictable way.) This is a kind of extremal dynamics, a scheme we've 
already encountered in the driven interface model; another familiar example is invasion per- 
colation ^^. Interestingly, the Bak-Sneppen model shows the same qualitative behavior as 



SODP: a singular stationary distribution of fitness values x*. The model exhibits avalanches 
in which replacement of a single species provokes a large number of extinctions. 

In the interface under extremal dynamics, the height Hi{t) cannot decrease. In the 
Bak-Sneppen model momentary setbacks are allowed {xj can decrease in a given step), but 
individuals of low fitness will eventually be culled. This is like an interface model with 
quenched noise such that, on advancing to a new position, an element may encounter a 
force that throws it backward, for a net negative displacement. The Bak-Sneppen model is 
equivalent to a driven interface in which the least-stable site and its neighbors are updated 
at the same moment; we can, as before, trade extremal dynamics for a limit of infinitely 
slow driving. 

Another way of obtaining the extremal dynamics of the Bak-Sneppen model as the 
limit of a stochastic process with purely local dynamics is as follows |^. Take a one- 



dimensional lattice (with periodic boundaries, for definiteness) , and assign random numbers 
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Xj, independent and uniform on [0,1], to each site j = 1, ..., L. The configuration evolves via 
a series of "fiips," which reset the variables at three consecutive sites. That is, when site j 
fiips, we replace Xj_i, Xj, and Xj^i with three independent random numbers again drawn 
uniformly from [0,1]. Let the rate of fiipping at site j be re"'^^^, where F"^ is a characteristic 
time, irrelevant to stationary properties. The Bak-Sneppen model is the (3 ^ oo limit of 
this process. 

We can get some insight into the stationary behavior via a simple analysis. Let p{x)dx 
be the probability that Xj G [x, x + dx]. The probability density satisfies 



dp{x) 
dt 



-e-^^'pix) - 2 / e-%(x, y)dy + 3 / e-^yp(y)dy (13) 

JO Jo 



where p{x, y) is the joint density for a pair of nearest-neighbor sites. If we invoke a mean-field 
factorization, p{x,y) = p{x)p{y), then 



dp{x) 
dt 

where 



The stationary solution is 



p{x) e-^" + 2/(/5) +3/(/3), (14) 



r(/3)^ f\-f^yp(y)dy . (15) 

J (J 



- 3 1 - e-'^f/^ 

Pstix) - 2 ^ _ ^_2^/3 ^ ^_^,^^^/3 _ ^^ . (16) 

The solution is uniform on [0,1] for /5 = 0, as we'd expect, but in the /3 ^ cxd limit we have 
Pst = (3/2)9(x — 1/3)9(1 — x). The probability density develops a step-function singularity, 
as in the Bak-Sneppen model. Not surprisingly, the mean-field approximation yields a rather 
poor prediction for the location of the singularity, which actually falls at 0.6670(1) |S^. (A 
two-site approximation places the singularity at x = 1/2.) The main point is that to realize 
singular behavior from a local dynamics, we have to tune a parameter associated with the 
rates. Alternative mean-field treatments of the Bak-Sneppen model may be found in Refs. 
5§ and H 



We can construct a model with the same local dynamics as that of Bak and Sneppen by 
replacing Xj_i, Xj, and Xj+i at rate 1, if and only if Xj < r. (Sites with Xj > r may only 
change if they have a nearest neighbor below the cutoff.) In other words, only sites with 
Xj < r are active; an updated site is active with probability r. There is an absorbing phase 



for small r, separated from an active phase by a critical point at some r^ [pO|-p^. To get 
the Bak-Sneppen model we forget about r, and declare the unique active site in the system 
to be the one with the smallest value of r. In the infinite-size limit, the probability to find a 
site with r < re is zero, in the stationary state. We see once again that in extremal dynamics 
we tune the order parameter itself to zero: at each instant there is exactly one active site, 
so pa = l/L. 

Grassberger and Zhang observed that the existence of SODP "casts doubt on the sig- 
nificance of self-organized as opposed to ordinary criticality." A similar doubt might be 
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prompted by our recipe for turning a conventional critical point self-organized. Of course, 
even if it is possible to explain all instances of SOC in terms of an underlying conventional 
critical point, the details of the critical behavior remain to be understood |]63|. Numerical 
results indicate that sandpiles, driven interfaces, and the Bak-Sneppen model define a series 
of new universality classes. Furthermore, no one has been able to derive the critical expo- 
nents of avalanches in SOC sandpiles, even the in abelian case, where quite a lot is known 
about the stationary properties [Q. 

V. SOC AND THE REAL WORLD 



Since SOC has been claimed to be the way "nature works" |^, we would expect to 
find a multitude of experimental examples where this concept is useful. Originally, SOC 
was considered an explanation of power laws, that it provided a means whereby a system 
could self-tune its parameters. So once we saw a power law we could claim that it was self- 
generated and "explained" by SOC. The previous sections should have convinced the reader 
that there are no self-tuning critical points, although sometimes the fine tuning is hidden, 
as in sandpile models. Therefore, an "explanation" of experimentally observed power laws 
requires the identification of the tuning parameters controlling the scaling, as in any other 
ordinary critical point. 

Here, we will restrict the discussion to experimental examples of avalanche behavior, 
leaving aside fractals and 1// noise whose connection with SOC is rather loose. (It is worth 
mentioning that a physical realization of self-organized criticality — without avalanches, 
as far as is known — has been identified in liquid ^He at the A point pH|.) Following the 



introduction of SOC, there were many experimental studies of avalanches, which sometimes 
yielded power-law distributions over a few decades, leading to endless discussions about 
the applicability of SOC. If we accept that self-tuned critical points don't exist, then these 
controversies have no basis: we have only to understand how far the system is from the 
critical point, and why. This task has only been accomplished in a few cases; several examples 
require further study, both experimental and theoretical. 

Soon after the sandpile model was introduced, several experimental groups measured the 
size-distribution of avalanches in granular materials. Unfortunately, real sandpiles do not 
seem to be behave as the SOC sandpile model. Experiments show large periodic avalanches 
separated by quiescent states with only limited activity |^|. While for small piles one could 



try to fit the avalanche distribution with a power law over a limited range |^S[, the behavior 
would eventually cross over, on increasing the system size, to the one described above, which 
is not scale- invariant. The reason sand does not behave like an ideal sandpile is the inertia 
of the rolling grains. As grains are added, the inclination of the pile increases until it reaches 
the angle of maximal stability 9s-, at which point grains start to flow. Due to inertia, the 
flow does not stop when the inclination falls to 6'^, but continues until the inclination attains 
the angle of repose 9s < 9c [^. Since the "constant force" (i.e., with 9 controlled) version 
of the system has a first-order transition, it is no wonder that criticality is not observed in 
the slowly driven case. So if we want to see power-law avalanches we have to get rid of the 
inertia of the grains. Grains with small inertia exist and can be bought in any grocery store: 
rice! A ricepile was carefully studied in Oslo: elongated grains poured at very small rate 
gave rise to a convincing power-law avalanche distribution 
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The previous discussion tells us that in order to observe a power-law avalanche distribu- 
tion, inertia should be negligible. As discussed in Sec. IV, the motion of domain walls in 
ferromagnets and flux lines in type II superconductors is overdamped, due to eddy-current 
dissipation; these systems are probably the cleanest experimental example of power-law dis- 
tributed avalanches. The noise produced by domain wall motion is known as the Barkhausen 
effect, first detected in 1919 fTlI]. Since then, it has become a common non-destructive 
method for testing magnetic materials, and its statistical properties have been studied in 
detail. When the external magnetic field is increased slowly, it is possible to observe well 
separated avalanches, whose size distribution is a power-law over more than three decades 
|7^^r6[. Domain walls are pushed through a disordered medium by the magnetic field, so 
we would expect a depinning transition at some critical field H = He- One should note, 
however, that the "internal field" acting on the domains is not the external field, but is 
corrected by the demagnetizing field H^ ~ —NM where M is the magnetization []75| , |76| and 
A^ the demagnetizing factor. Therefore, if we increase the external field at constant rate 
c, the internal field is given by Hint = ct — NM = ct — ky{t), where y{t) is the average 
position of the domain wall and k cc N. We recognize here the recipe for SOC given in 
section IVA: in the limit c —>■ and /c — > we expect to reach the critical point. This fact 
was indeed verified in experiments, where k can be controlled by modifying the aspect ratio 



of the sample [ 76 1 



In type II superconductors, when the external field is increased, flux lines are nucleated 
at the border of the sample and pushed inside by their mutual repulsion. The resulting 
flux density gradient, known as the Bean state |j7^, bears some analogy with sandpiles, as 
pointed out by De Gennes over 30 years ago [j78[. Unlike sand grains, flux lines have little 
inertia, and exhibit power-law distributed avalanches [79|. It is still unclear whether in this 



system a mechanism similar to the demagnetizing field maintains a stationary avalanche 
state, as in ferromagnets. Simulations of flux line motion ||80| have reproduced experimental 
results in part, but a complete quantitative explanation of the phenomenon is lacking. 

Another broad class of phenomena where SOC has been invoked on several occasions 

is that of mechanical instabilities: fracture, plasticity and dislocation dynamics. Materials 

subject to an external stress release acoustic signals that are often distributed as power laws 

over a limited range: examples are the fracturing of wood [^, cellular glass |8 

83| , in hydrogen precipitation |^ 



and concrete 
SHI. While it 



and in dislocation motion in ice crystals [|^ 
has often been claimed that these experiments provided a direct evidence of SOC, this is far 
from being established. In fact, fracture is an irreversible phenomenon and often the acoustic 



emission increases with the applied stress |8T|] with a sharp peak at the failure point. There 
is thus no stationary state in fracture, and it is debated whether the failure point can even 
be described as a critical point [Q or a first-order transition |8^. The situation might be 
different in plastic deformation, where a steady state is possible ||8^; recent experimental 
measurements of dislocation motion appear promising [Q. We may mention some related 
phenomena in which avalanches have been observed, and a theoretical interpretation is still 
debated: martensitic transformations [^, sliding systems |^ and sheared foams PT|] . 

Finally, it is worth mentioning that SOC has been claimed to apply to several other 
situations in geophysics, biology and economics. We have deliberately chosen to discuss 
only those examples for which experimental observations are accurate and reproducible. 
Even in these cases, it is often hard to distinguish between SOC-like behavior and other 



17 



mechanisms for generating power laws. This task appears almost hopeless in situations 



where only limited data sets are available, such as for forest fires [£^, or evolution ^3 
and remains very complicated in other cases, such as earthquakes, as witnessed by the vast 
theoretical literature on the subject [Q. 

VI. SUMMARY 

The genesis of self-organized criticality is a continuous absorbing-state phase transition. 
The dynamical system exhibiting the latter may be continuous or discrete, deterministic 
or stochastic, conservative or dissipative. To transform a conventional phase transition to 
SOC, we couple the local dynamics of the dynamical system to an external supervisor, or 
to a "drive" (sources and sinks with rates {h}). The relevant parameter(s) {(} associated 
with the phase transition are controlled by the supervisor or drive, in a way that does not 
make explicit reference to {(}. One such path involves slow driving {h —>■ 0), in which 
the interaction with the environment is contingent on the presence or absence of activity 
in the system (linked to {(} via the absorbing-state phase transition). Another, extremal 
dynamics, restricts activity to the least stable element in the system, thereby tuning the order 
parameter itself to zero. Specific realizations of this rather abstract (and general) scheme 
have been discussed in the preceding sections: sandpiles, forest fires, driven interfaces, and 
the Bak-Sneppen model. 

Viewed in this light, "self-organized criticality" refers neither to spontaneous or 
parameter-free criticality, nor to self-tuning. It becomes, rather, a useful concept for de- 
scribing systems that, in isolation, would manifest a phase transition between active and 
frozen regimes, and that are in fact driven slowly from outside. 
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Figure Captions 

Fig. 1. Stationary density p of active sites versus density of walkers C, in one-dimensional 
ARW. The inset is a logarithmic plot of the same data, where A = (^ — (^c- The slope of the 
straight line is 0.43. 

Fig. 2. Stationary avalanche-size distribution in the one-dimensional Manna sandpile with 
sequential dynamics, for L = 500, 1000, 2000, and 5000 (left to right) . Lower inset: finite- 
size scaling plot of the data in the main graph, InP* versus Ins*, with s* = L^^'^^s and 
P* = L'^'^^P. Upper inset: stationary density ( in the inner 10% of the system, plotted 
versus 1/L. The diamond on the ( axis is the critical density of ARW. 
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